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ABSTRACT 

We investigate the Luminosity - Metallicity {L — Z) relation in the local universe 
(0 < z < 0.15) using spectra of 6 387 star-forming galaxies extracted from the 2dF Galaxy 
Redshift Survey. This sample is by far the largest to date used to perform such a study. We 
distinguish star-forming galaxies from AGNs using "standard" diagnostic diagrams to build a 
homogeneous sample of starburst galaxies for the L — Z study. We propose new diagnostic di- 
agrams using "blue" emission lines ([O ll] Xilll, [O ill] A5007, and H/3) only to discriminate 
starbursts from AGNs in intermediate-redshift {z > 0.3) galaxies. Oxygen-to-hydrogen (0/H) 
abundance ratios are estimated using the "strong-line" method, which relates the strength of 
following bright emission lines [O ll] A3727, [O III] A5007, and H/3 (parameters R23 and O32) 
to O/H. We used the [N ll] A6584/Ha emission-line ratio as "secondary" abundance indicator 
to break the degeneracy between O/H and i?23. We confirm the existence of the luminosity - 
metallicity relation over a large range of abundances {^^ 2 dex) and luminosities (~ 9 magni- 
tudes). We find a linear relation between the gas-phase oxygen abundance and both the "raw" 
and extinction-corrected absolute B-band magnitude with a rms of ^ 0.27. A similar relation, 
with nearly the same scatter, is found in the R band. This relation is in good agreement with 
the one derived by Melbourne & Salzer (2002) using the KISS data. However, our L — Z re- 
lation is much steeper than previous determinations using samples of "normal" irregular and 
spiral galaxies. This difference seems to be primarily due to the choice of the galaxy sample 
used to investigate the L — Z relation rather than any systematic error affecting the O/H de- 
termination. We anticipate that this luminosity - metallicity relation will be used as the local 
"reference" for future studies of the evolution with cosmic time of fundamental galaxy scaling 
relations. 
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1 INTRODUCTION 

Our understanding of galaxy formation and evolution certainly ben- 
efits from improving our knowledge about the chemical properties 
of galaxies. The chemical composition of stars and gas within a 
galaxy depends on various physical processes, such as the star for- 
mation history, gas outflows and inflows, stellar initial mass func- 
tion, etc. Although it is a complicated task to disentangle the effects 
of these various processes, the determination of galactic chemical 
abundances at various epochs put strong constraints on the likely 
evolutionary histories of galaxies (see Pettini 2003 for a review). 
The correlation between galaxy metallicity and luminosity in 
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the local universe is one of the most significant observational re- 
sults in galaxy chemical evolution studies. Lequeux et al. (1979) 
first revealed that the oxygen abundance O/H increases with the 
total mass of irregular galaxies. To avoid several problems in the 
estimate of dynamical masses of galaxies, especially for irregulars, 
absolute magnitudes are commonly used instead. The luminosity - 
metallicity (L — Z) relation for irregulars was later confirmed by 
Skillman, Kennicutt & Hodge (1989), Richer & McCall (1995) and 
Pilyugin (2001a) among others. Subsequent studies have extended 
the L — Z relation to spiral galaxies (Garnett & Shields 1987; Zarit- 
sky et al. 1994; Garnett et al. 1997; Pilyugin & Ferrini 2000), and 
to elliptical galaxies (Brodie & Huchra 1991). The luminosity cor- 
relates with metallicity over ~ 10 magnitudes in luminosity and 2 
dex in metallicity, with indications suggesting that the relationship 
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may be environmental- (Vilchez 1995) and morphology- (Mateo 
1998) free. This suggests that similar phenomena govern the L — Z 
over the whole Hubble sequence, from irregular/spirals to ellipti- 
cals (e.g. Gamett 2002; Pilyugin, Vilchez & Contini 2004). Re- 
cently, Melbourne & Salzer (2002) have used a sample of 519 star- 
forming emission-line galaxies from the KPNO International Spec- 
troscopic Survey (KISS) to confirm the existence of the L — Z re- 
lation over a broad range of luminosity and metallicity. They foimd 
however that the slope of the L — Z relation is steeper than the 
dwarf galaxy L — Z relation (see also Pilyugin & Ferrini 2000). 
This may be evidence that the relationship is not linear over the full 
luminosity range (see also Mouhcine & Contini 2002). 

Many recent studies in galaxy evolution trace changes in scal- 
ing relations of galaxies as a function of cosmic epochs, such as 
the TuUy - Fisher relation for disks (e.g. Milvang-Jensen et al. 
2003; Ziegler et al. 2002) and the fvmdamental plane relation for 
spheroids (e.g. van Dokkum & Ellis 2003; Im et al. 2002). In this 
context, the luminosity - metallicity relation of galaxies can be 
used as a sensitive probe and consistency check of galaxy evolu- 
tion. 

The chemical properties of galaxies at different epochs pro- 
vide new constraints on theories of galaxy formation and evolution. 
If local effects such as the gravitational potential and "feedback" 
from supernova-driven winds are the dominant regulatory mecha- 
nisms for star formation and chemical enrichment, then the L — Z 
relation might be nearly independent of cosmic epoch, such as pre- 
dicted by the semi-analytic models of galaxy formation and evo- 
lution (e.g. Kauffmann ot al. 1993; Somcrvillc & Primack 1999). 
However, based on our current knowledge of galaxy evolution, sug- 
gesting that the cosmic star formation rate was higher in the past 
(e.g. Madau et al. 1996), and that the overall metallicity in the uni- 
verse at earlier times was lower, we might expect galaxies to be 
considerably brighter at a given metallicity if there was more pri- 
mordial gas available in the "young" universe to fuel star formation. 

With the advent of the 10-m class telescopes and the powerful 
optical and near-infrared spectrographs, it is now possible to probe 
the physical properties (SFR, extinction, chemical abundances, 
mass, stellar populations, etc) of intermediate (0 < z < 1; Kobul- 
nicky, Kennicutt & Pizagno 1999; Hammer et al. 2001, Contini et 
al. 2002, Lilly, CaroUo & Stockton 2003; Kobulnicky et al. 2003; 
Lamareille et al. in preparation) and high-redshift (1.5 < z < A; 
Pettini et al. 1998, 2001; Kobulnicky & Koo 2000; Mehlert et al. 
2002; Lemoine-Busserolle et al. 2003; Erb et al. 2003) galaxies. 
Even if the mmiber of galaxies per redshift bin is still small, these 
studies show that high-redshift galaxies {z > 2) are 2-4 magnitudes 
brighter than local galaxies of similar metallicity. This deviation 
from the local L — Z relation demonstrates that the ratio between 
the luminosity and metal content varies throughout a galaxy's life- 
time and is a powerful diagnostic of its evolutionary state. 

At intermediate redshifts, the situation is still not clear. Kob- 
ulnicky & Zaritsky (1999) and Lilly, Carollo & Stockton (2003) 
foimd their samples of intermediate- z galaxies to conform to the lo- 
cal L — Z relation without any significant evolution of this relation 
out to z ~ 1. On the contrary, Kobulnicky et al. (2003) claimed 
recently that both the slope and zero point of the L — Z relation 
evolve with redshift, the slope becoming steeper at early cosmic 
time. 

The determination of the fraction of galaxies that deviate from 
the local L — Z relation, as well as the amplitude of this deviation 
as a function of redshift, will certainly help in constraining the "dif- 
ferential" rate of star formation between high-redshift galaxies and 
local ones. However, such studies require a reliable local L — Z 



relation, extending over a large range of metallicities and luminosi- 
ties, taking into account the variety of the local galaxy population 
and giving a better estimate of the overall shape of the relation and 
its intrinsic scatter. 

The main goal of this paper is thus to establish the luminosity 
- metallicity relation in the local Universe. To do so we extracted 
a sub-sample of nearly 7 000 emission-line star-forming galaxies 
with high signal-to-noise spectra from the 2dF Galaxy Redshift 
Survey (hereafter 2dFGRS). We anticipate that this luminosity - 
metallicity relation will be used as the local calibration for future 
studies of the evolution with cosmic time of fvmdamental galaxy 
scaling relations. 

This paper is organized as follow. In Sect. 2, we present the 
selection criteria applied to the original 2dFGRS sample to define 
a sub-sample of emission-line galaxies with high quality spectra. 
Various observed line ratios and their implications for the analysis 
of the physical conditions in emission-line galaxies, and the oxygen 
abundance estimates are also discussed. In Sect. 3 we present the 
derived luminosity - metallicity relation, and show how it compares 
with previous determinations. In Sect. 4, we present the results of 
this paper and summarize our conclusions. 

Throughout this paper, all calculations assume the cosmol- 
ogy given by WMAP, with Qa = 0.73, Clm = 0.27 and Ho = 
71km s"^ Mpc"\ 



2 THE DATA 

2.1 Sample selection 

The sample of emission-line star-forming galaxies used to establish 
the local L — Z relation is extracted from the 2dFGRS observations. 

The 2dFGRS data set consists of optical (3600 - 8OOOA; 
AA = lOA) spectroscopy of more than 250 000 galaxies, covering 
two contiguous declination strips, plus 99 randomly located fields. 
One of the strips is located close to the super galactic plane, while 
the other strip is located in the northern galactic hemisphere. The 
mean redshift of the initial sample is 0.11, with almost all galaxies 
with z < 0.3. Full details of the survey strategy are given in Col- 
less et al. (2001). Although, the original 2dFGRS spectra are not 
flux caHbrated, the relative flux calibration, which is more critical 
for our analysis of emission-line ratios and equivalent widths is ac- 
curate (Lewis et al. 2002). The data reduction ensures that the flux 
ratio of adjacent emission lines is reliable. 

The measurement of equivalent widths has been done using a 
fully automatic procedure by fitting Gaussian profiles to both ab- 
sorption and emission lines simultaneously (see Lewis et al. 2002 
for a detailed discussion of the procedure and the determination of 
the fitting quality). The measured spectral features include [Oil] 
A3727, H/3, [Oiii] AA4959,5007, Hq, [Nil] A6584, and [Sii] 
A A67 17,6731. All strong emission lines needed to reliably estimate 
the nature of the main ionizing source and the gas phase oxygen 
abundance are measured. 

Starting from the original 2dFGRS sample of 269 013 spectra, 
we have applied different quality criteria to draw a sub-sample of 
7 402 galaxies. This sub-sample was selected as follows. 

To get accurate estimate of the gas phase abvmdances and to 
avoid any bias due to observational problems, we have restricted 
ourselves to galaxies with high-quality spectra. We first exclude all 
galaxies observed before 31th August 1999, which exhibit prob- 
lems due to a fault of the atmospheric dispersion compensator 
within the 2dF instrument (see Lewis et al. 2002 for more detailed 
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log([NII]A6583/Ha) log([SII]XX67 17,673 1/Ha) 

Figure 1. Diagnostic diagrams for our sub-sample of 7 353 narrow emission-line galaxies extracted from 2dFGRS. The continuous lines show the theoretical 
separation between starburst galaxies and AGNs from Kewley et al. (2001). An error of 0.15 dex has been added to the predicted separation in the [S II] 
AA6717-H6731/Ha diagram. The dashed line in the [N II] A6584/Ha vs. [O III] A5007/H/3 diagram shows the separation between starburst galaxies and AGNs 
as defined empirically by Kauffmann et al. (2003) using SDSS data. Dots represent star-forming galaxies, squares show AGNs, while triangles represent 
emission-line galaxies which are contradictory classified using both diagrams. 



discussion on this issue). For these galaxies, the fitting procedure 
to determine line properties leads to poor quality results. In addi- 
tion, only galaxies with good overall spectrum quality are selected. 
We thus selected galaxies for which we have the best quality fit, 
i.e. no bad pixels detected within 2 — 3cr of the line center, for 
emission lines used to determine both the nature of the main ion- 
izing source and the oxygen abundance, namely [O II] A3727, H/3, 
[Oiii] AA4959,5007, Ha, [Nil] A6584, and [Sii] AA6717,6731. 

We also choose spectra with a median signal-to-noise ratio 
(measured between 4000 — 7500 A) of at least 10. Indeed, we found 
during preliminary investigations that adding more noisy spectra 
will introduce a bias in the L ~ Z relation. Finally, for galaxies 
observed several times we kept the spectrum with the best overall 
quality or, if the same, the best signal-to-noise ratio. 



2.2 Nature of the main ionizing source: starburst or AGN? 

The selection criteria applied up to now focus only on the quality 
of the data. The restricted sub-sample built so far using these crite- 
ria includes emission-line galaxies with different ionizing sources, 
i.e., either young massive stars related to a recent starburst or a 
non-thermal continuum produced by an Active Galactic Nucleus 
(AGN). As we are interested in star-forming galaxies only, we have 
to exclude galaxies for which the emission-line spectrum is typ- 
ical of AGNs. AGN-like spectra are distinguished from starburst 
ones by stronger collisional emission lines relative to recombina- 
tion lines. 

Seyfert 1 galaxies are distinguished from the other types of 
AGNs (Seyfert 2, LINERs) by broader Balmer emission lines. The 
distribution of Full Width at Half-lVIaximum (FWHIVI) for [O II] 
A3727, [O III] A5007, and H/3 emission lines shows that the max- 
imum value for [Oil] A3727 and [Olll] A5007 is lOA, whereas 
49 galaxies have FWH]V[(H/3) exceeding this limit. These galaxies 



have been classified as Seyfert 1 galaxies, and thus excluded from 
the sample. 

The widely used technique to distinguish starburst galaxies 
from narrow-line AGNs invokes diagnostic line ratios (e.g. Bald- 
win et al. 1981; Veilleux & Osterbrock 1987). The large wavelength 
coverage of 2dFGRS spectra allows the easy identification of an 
AGN-dominated spectrum (Seyfert 2 and LINER) by the presence 
of high [Nil] A6584/HQ and [Sll] AA671 7-1-673 1/Hq flux ratios 
relative to [O III] A5007/H/3. 

Figure 1 shows the location of our galaxy sample in the di- 
agnostic diagrams [N II] A6584/Hq vs. [O III] A5007/H/3 and [S II] 
AA6717-H6731/Ha vs. [O III] A5007/H/?. We used equivalent width 
(EW) ratios instead of flux ratios in these diagrams. Note that these 
ratios are weakly sensitive to reddening and to uncertainties in the 
spectrophotometry. The left panel of Figure 1 shows that there are 
two well-separated sequences of emission-line galaxies in the [N II] 
A6584/HQ vs. [Olll] A5007/H/3 diagram. A similar distribution 
has been found for a sample of galaxies drawn from the Sloan Dig- 
ital Sky Survey (SDSS) and discussed recently by Kauffmann et 
al. (2003). Our sample contains galaxies with a large variety of 
excitation levels, suggesting that this sample contains both low- 
metallicity and metal-rich galaxies. This sample is thus suitable for 
studying the L — Z relation over a large range of metallicities. The 
AGN sequence originates from the bottom of the location of star- 
burst galaxies, i.e. where metal-rich galaxies locate. This suggests 
that the inclusion of AGNs in a sample used to investigate the L~Z 
relation will contaminate only the high-metallicity end. 

The location of starburt galaxies relative to AGNs in the di- 
agnostic diagrams has been extensively investigated in the litera- 
ture. Recently Kewley et al. (2001) have coupled stellar population 
and photo-ionization models to predict the distribution of starburst 
galaxies in the diagnostic diagrams. Theoretical model predictions 
have associated uncertainties. Several input parameters, such as 
massive star atmosphere models, stellar evolutionary tracks, deple- 
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tion factors and the slope of the initial mass function are not known 
with high accuracy. Observationally, starburst galaxies show cor- 
relations between their intrinsic properties, and consequently may 
scatter around any separation in the diagnostic diagrams. The solid 
line in the [S II] AA67 17+673 I/Ho vs. [Olll] A5007/H/3 diagram 
gives an indication of the model uncertainties, and represents an up- 
per limit to the theoretical boundary between starbursts and AGNs, 
corresponding to Kewley et al. (2001) model + 0.15 dex. 

In both diagnostic diagrams, the predicted separation of Kew- 
ley al. (2001) between star-forming galaxies and AGNs match quite 
well the observed separation between starbursts and AGNs. We 
also show, in the [N ii] A6584/Ha vs. [O ill] A5007/H/3 diagram, 
the separation limit adopted by Kauffmann et al. (2003) for SDSS 
galaxies. Unfortunately, they do not show either the location of 
SDSS galaxies in the [S ii] AA6717+6731/Ha vs. [O ill] A5007/H/3 
diagram or the separation between starbursts and AGNs in this di- 
agram. We found the separation limit adopted by Kauffmann et al. 
(2003) quite strict as it increases the number of galaxies with in- 
consistent classifications. We thus decided to adopt the theoretical 
seperation of Kewley et al. (2001). 

We classified galaxies as starburst or AGN only if the clas- 
sification criteria using both diagnostic diagrsims are consistent, 
and reject all galaxies with a different classification in each di- 
agnostic diagram. We end up with 7 085 starburst galaxies, 133 
narrow-line AGNs (113 Seyfert 2 and 20 LINERs) and 135 galax- 
ies with an ambiguous classification. The completeness of the final 
spectroscopicaly- selected sample is difficult to define. By restrict- 
ing the sample to the presence of all the emission lines, from [O ii] 
A3727 to [Sll] AA6717+6731, used both to classify the galaxies 
(starbiu'st/AGN) and to calculate the gas-phase oxygen abundances, 
we implicitly impose an upper redshift limit of 2: ~ 0.2. By also 
requiring high-quality data, we exclude all galaxies with redshifts 
higher than z > 0.15. Our sub-sample thus remains incomplete 
at higher redshift and fainter magnitudes. Since the analysis of the 
L — Z relation does not rely strongly on the completeness level of 
the galaxy sample, and in order to use the largest number of galax- 
ies for the analysis, we have decided to keep all the galaxies that 
satisfy our selection criteria listed above. Figure 2 shows the red- 
shift distribution of the 7 085 star-forming emission-line galaxies 
we will use to derive the local L — Z relation. The median redshift 
of this sample is < z >= 0.05. 

Taking advantage of the large sample of emission-line galaxies 
from 2dFGRS, we investigate the location of starbursts and narrow- 
line AGNs (Sey 2) in diagnostic diagrams involving "blue" emis- 
sion lines only, i.e. [O II] A3727, H/?, and [O III] A5007. Such di- 
agrams could be very useful to discriminate AGNs from starbursts 
in spectra of intermediate-redshift galaxies obtained in large ongo- 
ing or future deep spectroscopic surveys (VVDS, DEEP, etc). In- 
deed, the "red" emission lines, such as Ha, [N 11] A6584, and [S 11] 
AA67 17+6731, usually used to identify AGNs in the local universe 
are not observable in the optical range (up to ~ 9000 A) for red- 
shifts greater than ~ 0.3. Such investigations have already been 
performed in the past (e.g. Rola et al. 1997; Dessauges-Zavadsky 
et al. 2000) but on limited-size samples. 

The location of the "classified" narrow emission-line 2dFGRS 
galaxies in new diagnostic diagrams is shown in Figure 3. The first 
diagram between [Oiii] A5007/H/3 and [On] A3727/H/?, proposed 
by Rola et al. (1997), shows a clear distinction between starbursts 
(bottom left) and AGNs (upper right). From this diagram, we de- 
fine the following analytical expression for the demarcation curve 
between starbursts and AGNs: 
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Figure 2. Redshift distribution of 7 085 star-forming galaxies selected from 
the 2dFGRS sample. The mean value is 0.05. 
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Starburst galaxies are located below this line. The contami- 
nation by AGNs in this region is very low (~ 0.1%) taking into 
account ±0.15 dex uncertainties in the separation curve. 

The second diagram between R23 and O32 (see Sect. 2.4 for a 
definition of these ratios) shows also a clear separation between 
starbursts (left part) and AGNs (right part). From this diagram, 
we define the following analytical expression for the demarcation 
curve between starbursts and AGNs: 

log(032) = , ^ + 2.4 (2) 



log(R2 



1.7 



Starburst galaxies are located to the left of this line. The dis- 
crimation between starbursts and AGNs is even better using this di- 
agram, as the contamination by AGNs is very low (~ 0.1%) taking 
into account only ±0.10 dex uncertainties in the separation curve. 

We conclude that these two diagnostic diagrams could be very 
efficient to separate AGNs from starbursts in intermediate-redshift 
{z > 0.3) galaxies for which "blue" emission lines only are avail- 
able. 



2.3 Underlying Balmer absorption 

Balmer emission-lines (e.g. Ha, H/3) are generally affected by ab- 
sorption lines from the underlying stellar populations, older than 
the current starburst, and have to be corrected. Usually the effect of 
the underlying stellar H/3 in absorption is accounted for by apply- 
ing an ad hoc correction to the rest-frame equivalent width of the 
H/3 emission line independently of the galaxy properties. Different 
correction values for the stellar absorption have been used in the 
literature, ranging from lA (Jansen et al. 2000) to SA (Kennicutt 
1992). 

The H/3 emission lines in the spectra of our galaxy sample 
have been corrected from the underlying stellar absorption by a si- 
multaneous fit of the emission and absorption lines. Figure 4 shows 
the distribution of EW(H/3)abs restricted to galaxies with a good 
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log([0II]X3727/H^) log(R23) 

Figure 3. New diagnostic diagrams for our sub-sample of 7 353 naiTow emission-line galaxies extracted from 2dFGRS. The continuous lines show the proposed 
separation between starburst galaxies and AGNs (see text for analytical expressions). The dashed lines show the predicted separations plus an uncertainty of 
±0.15 dex and ±0.10 in the [0 III] A5007/H/3 vs. [0 II] A3727/H/3 and R23 vs. O32 diagrams respectively. Dots represent starburst galaxies, squares show 
AGNs, while triangles represent "unclassified" galaxies (see text for details). These new diagrams are very efficient to separate AGNs from starbursts in 
intermediate-redshift (2 > 0.3) galaxies for which "blue" emission lines only are available. 
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Figure 4. Distribution of the EW of H/3 in absorption restricted to galaxies 
with a good fit quality. The an'ow indicates the mean value SVF (H/3)abs = 
4.6A of the sample. 



quality fit to H/3 in absorption. The mean value of the EW(H/3)abs 
is 4.6A, similar to Kennicutt (1992) value. 

Instead of assuming a constant value for the EW of H/9 in 
absorption, as usually done in previous studies, we tried to link 
the EW(H/3)abs to easily measured galaxy parameters. We found a 
correlation with the (bj — rp) color, an indicator of galaxy morpho- 
logical type, i.e. redder is the galaxy smaller is the EW(H/3)abs- We 
found a lower dependancy on the absolute magnitude in bj-band 



(Mb ), i.e. brighter is the galaxy smaller is EW(H/3)abs- Figure 5 
shows the EW(H/3)abs versus (bj — rp) relation. Also shown is a 
linear least-squares fit to the relation described by the following 
empirical calibration: 



EW(H/3)abs = 7.33 - 0.908(bj - rp) + 0.0885]VIbj 



(3) 



The rms of the residuals is 1.75 A. The EW of Ha in ab- 
sorption is deduced from EW(H/3)abs using the empirical ratio 
EW(Ha)abB = 0.75EW(H/3)abs (Gonzalez-Delgado et al. 1999). 
For a number of galaxies in our sub-sample, i.e., 1 955, the fit for 
the H/3 absorption line is of a poor quality. These spectra have been 
corrected using Eq. 3. 



2.4 IVIetallicity 

Using estimates of the oxygen abundance as indicators of the gas- 
phase metallicity is now well documented and calibrated (e.g. Pagel 
1997). The most reliable method to derive the gas-phase oxygen 
abundance requires an estimate of the electronic temperature and 
density of the ionized gas (Osterbrock 1989). An accurate deter- 
mination of these parameters requires reliable measurements of 
temperature-sensitive auroral lines, usually the [O III] A4363 emis- 
sion line. Unfortunately, the 2dFGRS spectra do not have the re- 
quired S/N to measure correctly this line. The absence of direct de- 
tection of auroral lines, particularly in the case of metal-rich galax- 
ies where these lines are too weak to be observed, require the devel- 
opment of alternative methods based on strong emission lines. The 
most widely used method is based on the measurement of strong 
emission lines, such as [O II] A3727, [O III] AA4959,5007, and H/3. 
These lines contain enough information to get an accurate estimate 
of the oxygen abundance (McGaugh 1991). This is done through 
the parameter R23 introduced initially by Pagel et al. (1979), and 
defined as follow: 
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Figure 5. Variation of SW^(H/3)abs as a function of (bj — rp) color (top- 
panel) for our sample of starburst galaxies, restricted to galaxies with a good 
fit of H/3 absorption line. Solid, dotted and dashed lines show the empirical 
calibration found for Mb = —24, —18 and —12 respectively. Bottom- 
panel shows the residuals versus the (bj — rp) color, their rms is 1.75A. 

[OIII]AA4959, 5007 + [Oil] A3727 
R23 = ^ (4) 

Extensive studies have been dedicated to calibrate the relation 
between R23 and oxygen abundance (McCall et al. 1985; Pilyu- 
gin 2001b). The small number of emission lines needed to estimate 
R23 makes this method very attractive. Strong-line ratios reliably 
indicate the oxygen abundance to within the accuracy of the model 
calibrations, i.e. ±0.15 dex. Traditionaly, R23 is estimated from the 
flux of emission lines. Kobulnicky & Phillips (2003) have shown 
that the use of equivalent widths instead of fluxes of the same emis- 
sion lines to derive R23 gives almost the same results. The main 
advantage of this method is to be, at the first order, insensitive to 
the reddening. Thanks to the small wavelength separation between 
the emission lines involved in a given ratio, we are free to use the 
EWs of these lines instead of their fluxes. 

Generally, the [O in] A4959 emission line has a lower signal- 
to-noise ratio than the [Olll] A5007 line; thus the fitting proce- 
dure of this line gives poor results for ~ 40% of our galaxy 
sample. For these galaxies we have used the theoritical ratio 
[Oiii] A5007/[Oiii] A4959=2.85. The distribution of the [Oiii] 
A5007/[O III] A4959 ratio, restricted to galaxies with good fitting 
quality of both lines, is in good agreement with this value, taking 
into account measurement uncertainties. 

A complication with the use of R23 parameter to estimate the 
oxygen abundance is that the dependency of metallicity on this pa- 
rameter is degenerate. Indeed, at a fixed value of R23 two different 
values of metallicity are possible; different ionization parameters 
should lead to similar oxygen abundances. R23 increases with oxy- 
gen abundance in the low-metallicity regime (12 + log(0/H)< 
8.2), while for metal-rich objects (12 + log(0/H)> 8.4) it de- 
creases with O/H reflecting the efficiency of oxygen cooling over 
abundance in these objects. In the "intermediate" metallicity region 
(8.2 < 12 -I- log(0/H)< 8.4), galaxies may have a large range 
of metallicities for a tight range of R23. The uncertainties in this 



metallicity domain, i.e. whether an object with a given R23 param- 
eter lies on the metal-rich branch or on the metal-poor branch of the 
O/H vs. ^23 relation, dominate the tmcertainties related to model 
calibrations (see above). 

Different abundance indicators have been used to break the de- 
generacy, e.g. [Nil] A6584/[Oiii] A5007 (AUoin et al 1979), [Nil] 
A6584/[Oii] A3727 (McGaugh 1994), [Nil] A6584/Hq (van Zee 
et al. 1998), and galaxy luminosity (Kobulnicky, Kennicutt & Piza- 
gno 1999). The first indicator is sensitive to the ionization param- 
eter, while the galaxy luminosity can hardly be used to break the 
degeneracy, first of all because this is exactly what we are looking 
for, and secondly the universality of the L — Z relation is not es- 
tablished yet (see Contini et al. 2002 for detailed discussion). To 
break the degeneracy between low and high metallicities, we used 
the [Nil] A6584/Ha secondary indicator. Galaxies with log([Nll] 
A6584/Ha) < — 1 are classified as low metallicity objects, while 
the others are assumed to be metal-rich galaxies. 

The oxygen abundance can be determined using the calibra- 
tions of McGaugh (1991). Analytical expressions are found in Kob- 
ulnicky, Kennicutt & Pizagno (1999), both for the metal-poor and 
metal-rich regimes. These calibrations are parametrized as a func- 
tion of R23 and the ionization parameter defined as O32 = [O ill] 
AA4959,5007/[O 11] A3727. We found that 495 galaxies fall into 
the low-metallicity region and 5 892 into the high-metallicity one. 
For 698 galaxies, the oxygen abundance estimated using the low- 
metallicity caUbration was greater than the one derived using the 
high-metallicity calibration. This may happen when the R23 pa- 
rameter exceeds the maximum value allowed by the model. Those 
galaxies have been excluded from the sample in further analysis. 



3 THE LUMINOSITY - METALLICITY RELATION 

In this section, we investigate the luminosity - metallicity relation 
in the local universe (z < 0.15) using the sample of 6387 star- 
forming galaxies extracted from the 2dFGRS as explained in the 
previous section. 

Figure 6 shows our L — Z relation for star-forming emission- 
line galaxies. The general trend, widely discussed in the literature 
and confirmed here, is an increase of metallicty with luminosity 
over a large magnitude range, from M(bj ) — — 13 to —22. Both the 
least-square bisector fit (Isobe et al. 1990) and the average of the 
forward fit, considering the absolute magnitude as the independent 
variable, and the backward fit, considering the metalUcity as the 
independent variable, yield to similar results. The bisector linear 
fitt is shown as the soUd line in Figure 6, and is described by the 
following relation: 

12 + log(0/H) = 3.45(±0.09) - 0.274(±0.005) M(bj) (5) 

The rms scatter around the linear fit to the relationship be- 
tween metallicity and luminosity is equal to 0.27 dex , greater than 
the error estimates for metallicity determination (0.15 dex). The 
origin of this intrinsic scatter may be due to differences in the star 
formation history, the evolutionary status of the current starburst, 
different initial mass function, etc (see e.g. Mouhcine & Contini 
2002). The rms we found is identical to the one derived by Mel- 
bourne & Salzer (2002), who used a sample ^ 12 times smaller 
than our sample, suggesting that the rms value of 0.27 dex might 

t This fitting technique is more suitable in cases in which the intrinsic cor- 
relation being sought measurement errors are dominated by intrinsic scatter. 
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Figure 6. Luminosity - Metallicity relation for the sample of 6 387 star- 
forming galaxies extracted from the 2dFGRS. The linear regression on this 
sample is plotted as a solid line. The dotted line shows the relation from 
Melbourne & Salzer (2002), the short-dashed line the relation for normal 
galaxies (Kobulnicky, Kennicutt & Pizagno 1999), the long-dashed line for 
UV-selected and Hll galaxies (Contini et al. 2002) and the short-dashed- 
long-dashed line the relation for spirals from Pilyugin, Vilchez & Contini 
(2004). The bottom-panel shows the residuals of the linear regression. 

be the "real" scatter of the L ~ Z relation. The uncertainty on the 
coefficients of Equ. 5 is lower than the one derived by Melbourne 
& Salzer (2002) thanks to our larger sample. The linear regression 
estimated after eliminating the points distant for more than 3a from 
the original regression does not show any significant difference. 

It is instructive to compare our determination of the local 
L — Z relation with other published determinations. In Figure 6, 
we plot the L — Z relation for local "normal" irregular and spi- 
ral galaxies (Kobulnicky & Zaritsky 1999; short-dashed line), UV- 
selected and Hll galaxies (Contini et al. 2002; long-dashed line), 
spiral galaxies (Pilyugin et al. 2004, short/long-dashed line), and 
Melbourne & Salzer (2002) determination using a sample of 519 
emission-line galaxies from KISS project (dotted line). It is worth 
mentioning that all determinations of the L ~ Z relation have been 
corrected to the same cosmology. Our new luminosity - metallicity 
relation has a steeper slope than the one found for nearby dwarf 
irregulars (Skillman et al. 1989; Richer & McCall 1995). This is 
in contrast with Kobulnicky & Zaritsky (1999) and Garnett (2002) 
finding, who concluded that the L — Z relation exhibits a uniform 
growth over both the low- (irregular galaxies) and high-metallicity 
(spiral galaxies) regimes. On the other hand, our determination is 
in excellent agreement with the Melbourne & Salzer (2002) one. 
The slopes of both determinations are almost identical taking into 
account uncertainties. 

3.1 Possible sources of systematic errors 

3.1.1 Different methods to derive 0/H 

Before going further in our discussion, it will be useful to highlight 
an issue that may be problematic. A criticism that may be addressed 
is the suitability of the strong-line method to derive oxygen abun- 



dances, and to which extent the systematic errors inherent to this 
method may affect our L ~ Z relation. 

Different authors have pointed out that the R23, method in- 
volves systematic errors due to the failure to take into account 
the variety of physical conditions in different Hll regions. For in- 
stance, Kennicutt et al. (2003) suspected the _R23 method to be 
systematically biased toward higher value of 0/H than the direct 
method, especially at the high-metallicity end (see also Stasinska 
2002). The strong-line method yields oxygen abundances that are 
systematically higher than the electronic temperature-based tech- 
nique by an amount depending on the calibration and the exci- 
tation range considered. However, it is not clear whether the di- 
rect electronic temperature-based method is underestimating the 
oxygen abundances, or the strong-line method is overestimating 
them (see Stasinska et al. 2001 for a discussion of the modelling 
of Hll regions). The debate of understanding if the discrepancies 
in abundance scales are due to systematic biases in the electronic 
temperature-based scale or in the Hll region models used to cal- 
ibrate different strong-line vs. abundance relations, particularly at 
high metallicity regime, is not settled yet. 

The R2i parameter has been calibrated and correlated with 
metallicity by measuring oxygen abundances using the direct 
method, i.e. based on estimates of eletronic density and tempera- 
ture, for a large sample of Hll regions in the low-metallicity regime 
and using photoionization models for the high-metallicity regime 
(see Garnett 2002 for a review). The direct method suffers from 
severe limitations, such as the restricted range of the ionization pa- 
rameters and/or metallicities and the observability of intrinsic weak 
lines. This makes this method difficult to use for abundance es- 
timates, especially for galaxies at cosmological distances. There 
are several reasons for questioning the accuracy of both abun- 
dance scales (see Kennicutt et al. 2003 for a detailed discussion 
of different issues related to this topic). Pilyugin (2000) has pro- 
vided new calibrations to implement McGaugh (1991) corrections 
into the empirical strong-line method. This method, usually called 
the "P-method", is claimed to be more accurate for low-ionization 
Hll regions and for metal-poor nebulae than the traditional i?23 
method. However, the P-method seems to suffer some weaknesses 
even within the abundance range for which it was designed for (Sa- 
viane et al. 2002; Kennicutt et al. 2003). Edmunds & Paget (1984) 
and Pilyugin (2000) calibrations, between R2S, and 0/H, have the 
same slope but are shifted towards higher abundances by ~ 0.07 
dex, a value smaller than both the intrinsic uncertainties of the cali- 
bration and the intrinsic scatter of the luminosity - metallicity rela- 
tion. Moreover, Melbourne & Salzer (2002) have shown that using 
Pilyugin (2000) calibration instead of Edmunds & Paget (1984) one 
has only a small effect on the slope of the L ~ Z relation. 

To check any dependancy of the L — Z relation on the choice 
of the 0/H vs. J?23 calibration, we used the Pilyugin (2000) and the 
McGaugh (1991) calibrations to derive oxygen abundances. In each 
case, both line fluxes and equivalent widths are used to calculate the 
-R23 and 032 parameters. We have to mention that these calibra- 
tions are roughly parallel in the 0/H - i?23 diagram; the McGaugh 
(1991) calibration being systematicaly shifted toward higher abun- 
dances by ^ 0.2 dex at a given _R23- In each case, we performed 
a linear regression on the L ~ Z relation, in a similar way used 
to derive Equ. 5. We find that i) the slopes of the different L — Z 
relations are similar, ii) none of the adopted 0/H - _R23 calibration 
leads to a slope as shallow as the L — Z relation slope for irreg- 
ular galaxies. We conclude that while the slope of the luminosity 

- metallicity relation is slightly sensitive to the choice of the 0/H 

- -R23 calibration, our determination is not biased toward a steep 
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slope. As our first motivation is to build a luminosity - metallicity 
relation for the local universe to be used as a reference for inter- 
mediate and high-redshift studies where the majority of abundance 
estimates are based on the strong-line method (e.g. Lilly et al. 2003; 
Kobulnicky et al. 2003; Lemoine-Busserolle et al. 2003; Lamareille 
et al., in preparation), and since it is a "differential" comparison of 
the abundance properties of different galaxies, the accurate choice 
of the i?23 calibration is not a critical issue. 

3.1.2 Reddening 

Is there any other systematic effect that may be responsible for the 
steep L — Z relation we derived? The luminosity - metallicity rela- 
tion could be biased toward a steeper relation as a result of neglect- 
ing the effect of internal extinction. One may expect naively that 
bright galaxies are more affected by internal reddening than fainter 
ones, causing a systematic effect which increases as a function of 
the absolute magnitude and leading to a steeper L — Z relation. 
To investigate this effect, we have to estimate the extinction for our 
sample of galaxies. We derived the intrinsic interstellar reddening 
using the observed Balmer line flux ratios Ha/H/3 following the 
standard prescription (e.g. Osterbrock 1989), and using the extinc- 
tion law from Seaton (1979). Note that although different extinc- 
tion laws available in the literature are different in the UV, they 
show similar behaviour in the optical, making the results of the 
subsequent analysis independent of the chosen extinction law (dis- 
cussion of the galactic internal reddening and how its behaves as a 
function of galaxy properties will be discussed in more detail in a 
forthcoming paper). 

The results of the extinction correction on the L ~ Z rela- 
tion, using the O/H - i?23 calibration as given by Kobulnicky et al. 
(1999), is shown in Figure 7. Again, a linear least square bisector 
fit is applied to the data, and we find the following expression to 
the luminosity - metallicity relation: 

12 + log(0/H) = 4.07(±0.09) - 0.223(±0.004) M(bj) (6) 

with an rms of 0.32. The extinction-corrected L — Z rela- 
tion is slightly shallower than the previous one (uncorrected for 
extinction), but we still find a much steeper slope than what was 
reported before in the literature. Interestingly enough, the slope of 
our extinction-corrected L — Z relation is again in good agreement 
with the one found by Melbourne & Salzer (2002) after correcting 
their photometry for the internal extinction, despite the different 
methods adopted to estimate the internal reddening correction. The 
shift in the zero point between the two relations may be due to the 
different procedures used to correct for extinction. Melbourne & 
Salzer (2002) have corrected only galaxies that deviate at 2(t from 
the observed colour-magnitude relation, while we correct galaxies 
independently of their location in the observed colour-magnitude 
diagram. 

3.2 Discussion on the slope of the L ~ Z relation 

None of the possible sources of systematic errors listed above, i.e. 
the choice of the O/H - R23 calibration and the internal reddening, 
can explain the steep slope of the L — Z relation we find. How can 
we explain the difference between our relation and the one derived 
for "normal" irregulars and spirals? The fact that different deter- 
minations of the luminosity - metallicity relation for irregular and 
spiral galaxies, using different methods and techniques, have simi- 
lar slopes suggests that the systematic errors inherent to the O/H - 




M(bj) 

Figure 7. Extinction-corrected luminosity - metallicity relation for the sam- 
ple of 6 387 star-forming galaxies extracted from the 2dFGRS. The linear 
regression on this sample is plotted as a solid line. The dotted line shows 
the relation from Melbourne & Salzer (2002), the short-dashed line the re- 
lation for normal galaxies (Kobulnicky, Kennicutt & Pizagno 1999), the 
long-dashed line for UV-selected and Hll galaxies (Contini et al. 2002) and 
the short-dashed-long-dashed line the relation for spirals from Pilyugin, 
Vilchez & Contini (2004). The bottom-panel shows the residuals of the lin- 
ear regression. 

R23 calibration cannot bias the L — Z relation. Hence, the observed 
discrepancy between our relation and the other quoted determina- 
tions cannot be explained entirely by the uncertainties affecting the 
calibrations. 

The difference between these relations may be attribuable to 
the different samples used by various authors. We think that the 
nature of the galaxy sample used to study the luminosity - metal- 
licity relation has a strong impact on the derived slope. Mouhcine 
& Contini (2002) have investigated the L — Z relation for samples 
of Hll galaxies, UV-selected galaxies, and optical/far-infrared se- 
lected starburst nucleus galaxies. A linear fit to the entire sample 
provides a steeper correlation than what is found by Kobulnicky 
& Zaritsky (1999), i.e. a slope of 0.25. Once the starburst nucleus 
galaxies are excluded from the sample, the linear fit is consistent 
with shallower determinations of the L — Z relation (Richer & 
McCall 1985; Kobulnicky & Zaritsky 1999; Contini et al. 2002; 
Pilyugin et al. 2004). Again, this suggests that the L ~ Z rela- 
tion is primarily sensitive to the nature of galaxy samples. The 
linear fit to irregular and spiral galaxies shown in Figure 6 is al- 
most parallel to the fit to the UV-selected and Hll galaxies, i.e. at a 
given metallicity star-forming galaxies are brighter than "normal" 
galaxies. This is consistent with the expected effect of the current 
star formation episode in these galaxies, decreasing for a short pe- 
riod their mass-to-light ratio. As discussed by Melbourne & Salzer 
(2002), different recent determinations of the L — Z relation over 
a large range of metallicities/luminosities deviate from the relation 
derived for dwarf galaxies, the former being steepper (Zaritsky et 
al. 1994; Pilyugin & Ferrini 2000). The 2dFGRS galaxy sample 
contains a mixture of galaxies of different types, representative of 
the diversity of galaxies in the local universe. Thus, we believe that 
our determination of the luminosity - metallicity relation is con- 
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sistent and more representative of the L — Z relation in the local 
universe than the previous determinations restricted to a particular 
type of galaxies. 

The difference between our L — Z relation and the one de- 
rived for dwarf galaxies suggests that the overall form of the lu- 
minosity - metallicity relation may not be simply approximated 
by a single linear relation. We split our sample of galaxies into 
two sub-samples, the first one having 12 + log(0/H) < 8.3, as 
a representative sample of nearby metal-poor galaxies, and a 
second one having 12 -|- log(0/H) > 8.3 to mimic the metal- 
rich galaxy population. Figure 8 shows the L — Z relation 
for these sub-samples. Also shown are i) the linear fit to the 
metal-poor sample (solid line) given by the following equation: 
12 + log(0/H) = 4.29(±0.42) - 0.20(±0.02) M(bj) compared 
to Skillman et al. (1989) L ~ Z relation for nearby irregular 
(metal-poor) galaxies (dotted line), and ii) the linear fit to the 
metal-rich sample (dashed line) given by the following equa- 
tion: 12 + log(0/H) = 2.57(±0.32) - 0.30(±0.02) M(bj). The 
luminosity - metallicity relations derived for these sub-samples are 
different, with metal-rich galaxies following a steeper L — ^ re- 
lation than the metal-poor galaxies. The slope of the dwarf galaxy 
L — Z relation is slightly steeper than the one derived by Skill- 
man et al. (1989), for which the oxygen abundances were esti- 
mated using the direct electronic-temperature method. Restricting 
the 2dFGRS sub-sample to the best quality spectra, we found a 
slope (0.164 ± 0.02), for the dwarf galaxy L — Z relation, similar 
to what has beenderived by Skillman et al. (1989), i.e. a slope of 
0.153. The variation in the derived slopes being small and included 
in the error bars, it shows that the choice of the 2dFGRS sub-sample 
resulting from the selection criteria does not change significantly 
the slope of the derived L — Z relations. The small differences in the 
slope and zero point between our L ~ Z determination and the one 
by Skillman et al. (1989) for metal-poor galaxies may be due to: 
i) the different methods applied to derive O/H in low-mass galax- 
ies; Skillman et al. (1989) used the "direct" method as we used the 
strong-line one, //) the filter used to measure the galaxy magnitude 
varies also; 6j -filter in our case instead of the standard B-filter used 
by Skillman et al. (1989), and iii) the limited-size of the Skillman 
et al.'s sample compared with the 2dFGRS sub-sample. 



3.3 Dependance with other galaxy properties 

The large sample presented in this paper gives an opportunity to in- 
vestigate a possible sensitivity of the luminosity - metallicity rela- 
tion to various galaxy properties, such as star formation rate (SFR), 
galaxy stellar content, etc. 

At a given mass, starburst galaxies tend to be brighter than 
"normal" galaxies forming stars at a lower rate, especially in the 
B band. This could introduce a systematic shift in the L ~ Z re- 
lation for starburst galaxies toward higher luminosities for a given 
metallicity. Unfortunately, the 2dFGRS spectra are not calibrated 
in absolute flux, preventing us from estimating the current SFR of 
galaxies. Rather then using fluxes, we thus use Balmer emission- 
line EWs. The Balmer emission-line EWs are a measure of the rel- 
ative proportion of ionizing photons (produced by massive stars 
related to the current star formation event) and continuum photons 
(produced by a mix between the whole cluster embedded in the Hll 
region, the underlying older stellar population, and a contribution 
from the ionized gas). The Balmer emission-line EWs may be un- 
derstood as a measurement of the specific SFR. An analysis of the 
luminosity - metallicity relation as a function of the EW of Balmer 




Figure 8. Luminosity-metallicity 

relations for 2dFGRS metal-poor (12 + log(0/H) < 8.3; tiiangles) and 
metal-rich (12 -|- log(0/H) > 8.3; dots) galaxies. The solid line shows 
a linear fit to metal-poor galaxies, dashed line shows the fit for metal-rich 
galaxies, and the dotted line shows the L — Z relation of dwarf irregulai' 
galaxies by Skillman et al. (1989). The bottom-panel shows the residuals of 
the linear regression. 

emission lines might therefore tell us something about the sensitiv- 
ity of this relation to the "current" SFR of galaxies. 

We split the sample into different groups as a function of their 
EW(Ha)e. The 
first group contains galaxies having EW(HQ:)e > 40A, group II 
contains galaxies with 20A < EW(Ha)o < 40A, group III con- 
tains those with lOA < EW(Hq)c < 20A, and finally group IV 
contains galaxies having EW(Hq)c < lOA. We then perform, for 
each group, a linear least square bisector fit to the L — Z relation. 
We find that the luminosity - metallicity relations for the first two 
groups are similar, whithin the errors, to each other and to the whole 
sample L ~ Z relation, showing similar scatter around the fit. This 
might be expected as the first two bins are sampling a short starburst 
age interval (Stasinska et al. 2001), within which the luminosity is 
not expected to vary significantly. The two remaining groups show 
L~ Z relations with different slopes, indicating a sensitivity of the 
latter to the current SFR. However, the scatter around the fits are 
large, preventing us from drawing any firm conclusion. 

Similarly, we investigate the sensitivy of the luminosity - 
metallicity relation to the galaxy colour. Again, no firm conclusion 
can be drawn. Furthermore there is no obvious correlation between 
the residuals from the L — Z relation in one hand and EW(Ha)c 
and colour in other hand. The present sample of galaxies does not 
present any appealing evidence for the sensitivity of the luminosity 
- metallicity relation to the stellar content of galaxies. 

3.4 The L ~ Z relation in the R band 

rp-band photometry is available for the 2dFGRS galaxy sample. 
We thus investigate the rp-band luminosity - metallicity relation 
and find a trend between these two quantities similar to what is 
found in the "blue" band. We fit a linear relation between the oxy- 
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Figure 9. Luminosity-Metallicity relation for different bins of redshift: z < 
0.05 (upper-left, blue dots), 0.05 < z < 0.1 (upper-right, green crosses) 
and z > 0.1 (bottom-right, red triangles). The bottom-left panel show the 
three redshift bins together The solid line shows the L — Z relation given 
in Eq.5. 



gen abundance and the rp-band absolute magnitude, in a similar 
way to what has been done for the bj-band relation. We find the 
following rp-band L ~ Z relation: 



12 + log(0/H) =3.72(±0.l 



0.249(±0.004) M(rp) (7) 



with an rms of the residuals of 0.25 dex, lower than the scatter 
for the 6j-band L — Z relation. This might be due to the smaller 
sensitivity of rp-band luminosity to the occurrence and the evolu- 
tionary status of a recent star formation episode. 



3.5 Dependance with redshift and selection biases 

In this last section, we investigate whether the L — Z relation is 
redshift-dependant as recently suggested in Schulte-Ladbeck et al. 
(2003) using SDSS data. Despite the small redshift range of our 
sub-samble (0 < z < 0.15), Figure 9 shows indeed significant dif- 
ferences if we divide the sample into three different redshift bins 
(i.e. z < 0.05, 0.05 < z < 0.1 and 2: > 0.1). At a given metal- 
licity, the available observed luminosity range seems to be shifted 
toward higher luminosities when we look further in the past. This 
effect is supported by calculating, for each bin, the mean value of 
the luminosity, which increases (—18.2 for z < 0.05, —19.2 for 
0.05 < 2: < 0.1 and —20.0 for z > 0.1) while, in average, the 
metallicity remains approximatively constant (12 + log(0/H) = 
8.58, 8.70 and 8.77 respectively). This effect is particulary strong 
for low-metallicity objects (12 -I- log(0 /H) < 8.3, see the bottom- 
left panel). 

Before drawing any conclusion on the evolution of the L — Z 
relation with redshift, one must take into account some selection 
effects. First of all, we have to consider that high-luminosity ob- 
jects will preferably appear at high redshifts, where the volume of 
the observed universe is large enough to make their number signi- 
ficative given their small intrinsic proportion. In the opposite way, 
low-luminosity objects may appear at all redshifts, but most of them 



are observed at low redshift due to the magnitude-limited nature of 
the survey (the Malmquist bias). The first effect does not explain 
the evolution of the mean luminosity, while it is normalized to the 
total number of observed objects. The last one is more problematic 
because we loose almost all low-metallicity objects at higher red- 
shift. Our quality selection criteria also introduce a bias: very metal 
rich, and thus very luminous, objects have been rejected from our 
sub-sample by the quality criterium on oxygen emission lines, con- 
sidering that they are very faint (low S/N) for 12 + log(0 /H) > 9. 

Taking into account these selection effects, no firm conclusion 
can be drawn from this study on any evolution of the L — Z relation 
with redshift. 



4 CONCLUSIONS 

We have investigated the luminosity ~ metallicity relation in the 
local universe, using data of ~ 7000 star-forming galaxies (0 < 
z < 0.15) extracted from the 2dFGRS spectroscopic data set. The 
sample used in this paper is by far the largest sample of galaxies to 
date to derive the L — Z relation. 

We first distinguished star-forming galaxies from AGNs us- 
ing "standard" diagnostic diagrams to build a homogeneous sam- 
ple of 6 387 starburst galaxies for the L — Z study. We propose also 
new diagnostic diagrams using "blue" emission lines ([O II] Xilll, 
[O III] A5007, and H/3) to discriminate starbursts from AGNs in 
intermediate-redshift {z > 0.3) galaxies. We have shown that the 
underlying absorption lines affecting Balmer lines in emission de- 
pend on galaxy properties, with a wide distribution ranging from 
1 A to 10 A with an average value of 4.6A. Oxygen-to-hydrogen 
abundance ratios were estimated for this sample of star-forming 
galaxies using the strong-line method, which relates the R23. and 
032 parameters to 0/H. We confirmed the existence of the lumi- 
nosity - metallicity relation over a large range of abundances and 
luminosities. We find a linear relation between the gas-phase oxy- 
gen abundance and both the "raw" and extinction-corrected 6j-band 
absolute magnitude with a rms of ~ 0.26. A similar relation, with 
nearly the same scatter, exists for the rp-band absolute magnitude. 

Our determination of the L — Z relation is in a remark- 
ably good agreement with the one derived by Melbourne & Salzer 
(2002) using the KISS data . However, the slope of the relation we 
derived is much steeper than the ones previously determined using 
different samples of "normal" dwarf and spiral galaxies. We argue 
that this difference is not due to any systematic error inherent to the 
method used to derive O/H in galaxies. The nature of the galaxy 
sample used to investigate the L — Z correlation is however cru- 
cial. Neglecting a certain type of galaxies, such as starburst nucleus 
galaxies, in previous studies has biased the determination of the 
L — Z relation toward shallower slopes. The L — Z relation re- 
stricted to metal-poor galaxies of our sample is in agreement with 
the relation derived by Skillman et al. (1989) for the low-metallicity 
irregular galaxies. 

The luminosity - metallicity relation extends over a range of 9 
magnitudes in luminosity, and a factor of ~ 50 in oxygen abun- 
dance. The 2dFGRS sample of galaxies used in this paper con- 
tains a large diversity of masses and stellar populations, from dwarf 
galaxies to massive spirals. By using such a sample of galaxies, 
we are constructing a more general luminosity - metallicity rela- 
tion in the local universe than previous studies restricted to a given 
class of galaxies. This correlation seems to be continuous from 
faint dwarf galaxies to massive spirals, implying that the phys- 
ical mechanism(s) regulating this correlation is common to star- 
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forming galaxies over the whole Hubble sequence, without a clear 
dependency on the stellar content of galaxies. 

The investigation of the chemical properties of high-a galax- 
ies, and the comparison between galaxy properties over a large 
range of redshifts will certainly help to clarify the evolutionary pat- 
tern of galaxies at different cosmic epochs. 
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